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Appendix  A 
DATA  REPORT 

This  appendix  summarizes  the  wake  vortex  and  wind  measurements 
made  with  the  Lockheed -Huntsville  LDV  during  the  research  program.  The 
logs  for  the  LDV  measurements  showing  the  run  identification,  time  of  day, 
estimated  wind,  and  operator  comments  are  presented  on  pages  A -2  through 
A- 11.  The  cataloged  wake  vortex  measurements  are  presented  on  pages 
A- 12  through  A-253  in  terms  of  the  trajectory  of  the  wake  vortex  with  re- 
spect to  the  runway  centerline.  At  the  start  of  each  new  tape,  the  input 
parameters  selected  for  the  VAD  and  Vortex  Track  program  are  ven 
for  the  wake  trajectory  calculations. 
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DATA  LOG  OF  WAKE  VORTEX  AND  WIND  VELOCITY  MEASUREMENTS 
CONDUCTED  BY  LOCKHEED-HUNTSVILLE  AT  JFK 


Tape 

ID 

JFK 

Run  ID 

Time 

Est.  Wind 
Azimuth 
(from) 

Comments  j 

No. 

AC  Type 
or  VAD 

Start 

Stop 

! 

1 

1 

VAD 

10:51 

10:57 

NW 

120°  total  inc.  angle 

1 

2 

VAD 

10:03 

11:06 

NW 

60°  total  inc.  angle 

2 

1 

VAD 

03:50 

03:54 

NW 

Very  clear  dav;  VAD  scan  heights 

50.  60,  70,  80.  100,  200,  250 

2 

2 

747 

03:56 

03:57 

NW 

Changed  at  three 

? 

3 

707 

03:58 

04:00 

NW 

4 

727 

04:01 

04:02 

NW 

1 

i 

5 

DCIO 

04:08 

04:09 

NW 

20  dB  input  atten. 

6 

DC9 

04:12 

04:13 

NW 

1 

2 

7 

707 

04:16 

04:17 

NW 

2 

8 

727 

04:18 

04:19 

NW 

1 

O 

9 

707 

04:20 

04:21 

NW 

i 

2 

10 

707 

04:23 

04:24 

NW 

2 

1 1 

DCIO 

04:24 

04:25 

NW 

2 

1 2 

747 

04:34 

04:35 

NW 

Wind  at  or  near  run  end;  Teflon  off 

2 

13 

VAD 

NW 

2 

14 

DC8 

04:37 

04:38 

NW 

y 

15 

DCIO 

04:39 

04:39 

NW 

16 

DC9 

04:41 

04:42 

NW 

2 

1 7 

747 

04:43 

04:44 

NW 

2 

18 

747 

04:45 

04:46 

NW 

2 

19 

DC8 

04:47 

04:47 

NW 

1 

1 2 

20 

707 

04:49 

04:50 

NW 

1 

21 

VAD 

04:55 

05:00 

NW 

See  run  1 

L..  ' 

22 

VAD 

05:00 

05:03 

NW 

Parity  in  21st  record 

s 


A -2 


Run  ID  Time 


AC  Type]  Start  Stop 
of  VAD  I 


15:00:26  1 5:05:40 


747 
707 
74  7 
747 
747 
737 
747 
747 
727 
707 
707 
707 
747 
747 
747 
DC-9 
VAD 


DC-9 

707 

707 

DC-8 

707 

727 

727 

707 


Est.  Wind 
Azimuth 
(from) 

Comments 

15:09:35 
15:11:50 
15:13:40 
1 5:24:00 
I 5:28:00 
15:3  5:00 
15:37:05 
i 1 5:40:30 
1 5;44;40 
15:46:00 
1 5:48:55 
15:51:00 
j 15:52:15 
j 15:55:15 
I 15:58:30 

I 16:01:00 

I 16:04:45 


16:1  1:30 
16: 14:30 
16:17:45 
16: 19:35 
16:23:30 
16:26:00 


727  , 16:28:00 

707  I 16:30:00 

DC-  10  16:33:00 


1 5:10:35 
15:12:15 
15: 14:45 
15:24:55 
15:24:55 
15:35:45 
15:38:05 
1 5:41:35 
15:45:30 
15:47:00 
15:49:40 
15:51:40 
15:53:10 
15:56:30 
15:59:30 
16: 14:00 
16:09:00 


16:12:10 

16:16:15 

16:18:40 

16:20:45 

16:25:00 

16:27:15 

16:29:00 

16:31:10 

16:34:20 


rvw 

NW 

NW 

NW 

NW 

NW 

NW 

NW 

N-NE 

N-NE 

N-NE 

N-NE 

N-NE 

N-NE 

N-NE 

N-NE 

N-NE 


N-NE 

N-NE 

N-NE 

N-NE 

N-NE 

N-NE 

N-NE 

Calm 

Calm 


40,  45,  60,  75,  100,  150,  200. 
300  m at  1 rev.  per  altitude 


1 2 knots  at  340 
Looks  good 


360  at  knots 


< 5 knots 
Calm 


Good!  Low!  Calm  winds 


See  Run  1 notes  for  steps 
observed  wind  at  ^ 3 mph 
horiz.  comp.  707  passed 
16:04:07  about  record  57 

About  Run  1 9 

Good  run ! 

Freq.  going  about  F.  S. 

360°  at  8.  Good  one 
Good  one 
Good  one 

Good  one 


AC  Type 
of  VAD 


Stop 


Est.  Wind 
Azimuth 
(from) 


Comments 


Start 


727 

16:35:00 

16:36:00 

Calm 

Good  one,  came  in  high 

DC- 10 

16:37:00 

16:38:30 

Calm 

Good 

707 

16:40:00 

16:41: 10 

Calm 

Good 

747 

16:44:00 

16:45:00 

Calm 

Good 

DC-10 

16:46:30 

16:47:30 

Calm 

Good 

DC-  10 

16:49:00 

16:50:10 

Calm 

747 

16:51:30 

16:53:30 

Calm 

Calm.  Good!  Real  Good 

747 

16:54:00 

16:55:00 

Calm 

VAD 

17:57:00 

17:00:00 

Calm 

See  run  1 notes  for  steps 

L-10-1 1 

17:02:00 

17:03:30 

Calm 

L-10-1 1 

17:05:00 

17:07:00 

Calm 

Calm 

707 

17:11:15 

17:12:10 

Calm 

Calm 

727 

17:13:15 

17:14:00 

Calm 

Abort  run 

747 

17:15:00 

17:16:30 

1 

Calm 

Exceeded  max.  freq. 

707 

17:18:00 

17:19:00 

DC-8 

17:25:00 

17:26:20 

0080  at  5 knots 

747 

17:27:00 

17:29:00 

Real  good;  hanging  high 

727 

Mis  se( 

1 TM 

707 

17:32:00 

17:33: 30 

727 

17:35:00 

17:35:30 

007  at  5 

707 

1 7:37:00 

17:38:15 

007  at  5 

Wheel 

10:04:00 

10:09:00 

222  m for  1st  30  or  31 
records;  322  m for  rest. 
Wheel  range  (ft)  to  side  of 
van  (542) 

W heel 

1 0:33:00 

1 

1 

10:38:00 

1 

Wheel  at  328°  to  side  of 
van.  A/C  turned  on  off  at 
record  29  and  on 

Wheel 

1 hl5:00  1 

1 1;20:00| 

230°  distance  from  wheel 
to  van  side 

Abort 
W hcf] 


1 1:52:00 
1 1:55:00 


1 1:63:00 
12:00:00 


30  kHz  BW  to  RE  35 


Comments 

(1)  Scanner  running  clockwise 
looking  down  on  van  (JFK 
ID  5);  (2)  Conditions  clear; 
gusty  winds  toward  van  fror 
W/SW  (JFK  ID  6) 


Calib  1 12:02:00  I 

Calib  2 
Calib  3 

Calib  4 12:03:00 

VAD  14:10:00  14:15:00 

14:15:00  14:25:00 

I 14:25.00  14:40:00 

I 14:41:00  ! 14:46:00 
1 1 

I 14:46:00  I 14:56:00 

|14:56:00  15:1  1:00 

I 15: 12:  00  15:17:00 

I 15:17:00  i 15:27:00 
I I 

I 1 5:27:00  i 15:42:00 


230  at  10 
Using  226 
230  at  6 


200  at  5 
220°  at  5 
220  at  3 


VAD 

15:  12:00 

15:17:00 

SW 

747 

1 5:34:00 

1 5.3  5:00 

SW 

707 

1 5:40:00 

15:41:00 

290°  at  16 

707 

1 5:43:00 

15:44:00 

290°  at  16 

727 

1 5:46:00 

15:47:00 

290°  at  14 

707 

1 5:48:00 

15:49:00 

DC -9 

15:50:00 

15:51:00 

747 

1 5:52:00 

15:53:00 

290  at  16 

747 

1 5:54:00 

15:55:00 

707 

1 5:55:00 

15:56:00 

290  at  15 

35.  45,  60,  75.  100  150,  200, 
400  m at  1 rev/alt 

2 rev/alt 


4 rev/alt 
4 rev/alt 
2 rev/alt 
1 rev/alt 

1 rev/alt 

2 rev/alt 
4 rev/alt 
4 rev/alt 
2 rev/alt 
1 rev/alt 

Alt:  35,  45,  60,  75,  100.  1 50, 
200,  400  m at  1 rev/alt 

(280°  at  18)  azimuth  200°; 
data  bad,  need  280? 


Good  run 


Vortices  blowing 
in  toward  van 


i^st.  Wind 
Azimuth 

stop  i 


Comments 


11:01:00  I 
11:02:00  j 
1 1:03:00  ■ 
1 1:04:00  I 
1 1:05:00  | 

1 1:06:00  i 

1 1:07:00  | 

1 1:08:00  I 

1 1:09:00  I 

1 1:10:00  I 

11: 1 1:00  I 
1 1:12:00 
1 1; 14:00 
1 1: 1 5:00 
1 1: 16:00 
1 1: 17:00 
1 1:18:00 
1 1:19:00 


11:01:20 
1 1 :02:20 
1 1:03:20 
1 1:04:20 
1 1:05:20 
1 1:06:20 
1 1:07:20 
1 1:08:20 
1 1:09:20 
1 1: 1 0:20 
1 1:1  1:20 
1 1: 12:20 
11:14:20 
1 1: 1 5:20 
1 1:16:20 
1 1:17:20 
1 1:18:20 
1 1: 19:20 


1 1:20:00 
I 1:21:00 
1 1:22:00 
1 1:23:00 
I 1:24:00 
1 1:25:00 
1 1:26:00 
1 1:27:00 
1 1 :28:00 
1 1 ;29:00 


1 1:20:20 
1 1:21:20 
1 1:22:20 
1 1:23:20 
1 1:24:20 
1 1:25:20 
1 1:26:20 
1 1:27:20 
1 1:28:20 
1 1:29:20 


Ranui-  Scan 


Tape 

ID 

Run  ID 

T imo 

Est.  Wind 

Comments 

JFK 

No. 

AC  Type 
of  VAD 

Start 

Stop 

(from) 

8 

37 

RC" 

1 1:30;  00 

1 1:30:20 

8 

38 

1 1:31:00 

1 1:31:20 

8 

39 

1 1:42:00 

1 1:45:00 

8 

40 

1 1:49:30 

1 1:48:30 

8 

41 

VAD 

1 1:91:00 

1 1:52:30 

200°  at  7 

Halt  after  1}  min,  A^inds 

> full  scale 

8 

4 2 

VAD 

* 

RC 

1 1 :95:00 

1 2:00:00 

200°  at  7 

Shear  running  about  1.2 
m cycles;  gw  about  4U0  kHz 

8 

43 

12:08:30 

12:17:00 

Azimuth,  310°;  left  of  Tower  2 

8 

44 

12:44:00 

12:47:00 

Azimuth,  240°;  to  right  and  top. 

T ow  e r 1 

8 

45 

VAD 

12:56:00 

13:01:00 

r- 

0 

o 

o 

9 

VAD 

16:29:30 

16:31:30 

230  at  1 1 

Alt:  35,  41 , 50,  60,  70,  100, 
200,  300  m at  1 rev/alt 

9 

2 

DC -9 

16:41:00 

16:42:00 

230  at  1 1 

Gone  in  a hurry 

Azimuth  = 

9 

3 

DC  - 1 0 

1 6:44:00 

16:45:00 

Gone  in  a hurry 

280° 

9 

4 

707 

16:47:00 

16:49:00 

Azimuth  = 100° 

Notice  azi- 
muth change 

9 

747 

16:90:00 

16:51:00 

Good 

Good  display  s , 

9 

6 

DC- 10 

16:52:00 

16:53:00 

Good 

but  very  short 

9 

707 

16:54:00 

16:55:00 

Good 

lived. 

9 

DC- 10 

16:55:00 

16:56:00 

Good 

1 I 

1 

VAD 

1 1:43:00 

1 1:55:00 

220° 

Alt:  41  , 50,  60.  70.  100.  200. 

300  (SR  1) 

1 1 

) 

c. 

VAD 

1 2:46:00 

12:59:00 

200 

Ditto 

1 1 

3 

1 3: 1 5:00 

Ditto  plane  crash  at  1310 

1 1 

4 

1 3:59:00 

14:45:00 

300  m - 1.4  MH2 

1 I 

5 

14:52:45 

350  - 600  MH2 

70  - 600  MI-U 
1 00 . two  poaks 
2010  up  to  1 .8  Mllz  at  200” 


rape 

ID 

Run  ID 

T ime 

IFK 

No. 

— 

AC  type 
of  VAD 

Sta  rt 

Stop 

13 

1 

V.A.D 

14:50:00 

14:52:00 

1 5 

VAD 

14:55:00 

15:60:00 

VAD 

15:20:00 

15:25:00 

1 3 

4 

VAD 

1 5:48:00 

15:54:00 

1 3 

3 

VAD 

16:19:00 

16:24:00 

1 5 

6 

VAD 

16:49:00 

16:55:00 

1 3 

7 

VAD 

10:41:00 

10:52:00 

13 

8 

707 

1 1:27:00 

1 1:28:00 

1 3 

9 

DC-9 

1 1:34:00 

1 1:36:00 

13 

1 0 

300 

1 1:46:00 

1 1:48:00 

1 3 

n 

VAD 

1 1:52:00 

1 1:55:00 

1 3 

1 2 

DC-8 

1 1:58:00 

1 1:10:00 

1 3 

1 3 

707 

1 2:02:00 

12:04:00 

1 3 

14 

727 

12:08:00 

1 2:09:00 

1 3 

1 5 

727 

1 2: 1 0:00 

12:1  1:00 

1 3 

16 

DC -9 

12: 13:00 

1 2: 14:00 

1 3 

17 

707 

12:16:00 

12:17:00 

14 

1 

VAD 

14:32:00 

14:34:00 

z 

747 

14:39:00 

14:40:00 

14 

3 

VAD 

14:49:00 

14:52:00 

14 

707 

14:54:00 

14:55:00 

14 

5 

747 

14:56:00 

14:58:00 

14 

6 

747 

14:59:00 

15:00:00 

14 

7 

727 

1 5:01:00 

15:02:00 

14 

8 

DC -9 

15:03:00 

1 5:  14:00 

14 

9 

L-10-11 

1 5:06:00 

1 5:08:00 

14 

10  ! 747 

1 5;  10:00 


15:12:00 

Est.  Wind 
Azimuth 
(from) 


Comments 


220 


220 

220'^ 

220 

220'^ 

220 

330*^ 


Alt;  35.  41,  50.  60,  70,  100, 
200,  300  m i 


o 


Alt:  35,  41 . 50.  60.  70,  100. 
200,  300 


Airbus  (France) 


Alt:  35,  41,  50.  60,  70.  100 
200,  300 


370  at 
330  at 


10 

13 


15 


1 1 


No  gas  — change  prg.  ! 

Alt:  35,  41  , 50.  60.  70.  100. | 
200,  300  m at  1 rev/alt  I 


Good 


Looks  pretty  good 
Looks  good 


'I  apo 
ID 

Run  ID 

me 

Est.  Wind 

Comments 

JFK 

No. 

.4.C  Type 
of  VAD 

Start 

Stop 

14 

1 1 

707 

15:14:00 

15:15:00 

O.  K. 

14 

1 2 

DC -8 

1 5:  16:00 

15: 17:00 

1 

14 

1 3 

707 

15:20:00 

15:21:00 

Abort,  see  nothing,  (work- 
ing on  processor)  | 

14 

14 

707 

1 5:22:00 

15:22:00 

Abort;  nothing 

14 

15 

707 

15:25:00 

15:26:00 

Notice  30  kHz  bandwidth 

14 

16 

707 

1 5:34:00 

1 5:35:00 

14 

17 

707 

15:35:00 

15:36 

Looks  good 

14 

18 

707 

15:36:00 

15:37:00 

330  at  10 

O.  K.  ! 

14 

19 

7 37 

1 5:38:00 

1 5:39:00 

Best  yet  j 

14 

ao 

7 07 

15:40:00 

15:41:00 

340  at  1 2 

Good  j 

14 

ji 

707 

15:42:00 

15:43:00 

1 

14 

22 

DC-9 

15:44:00 

1 5:44:00 

Not  good  , 

.4 

23 

DC-9 

15:46:00 

1 5:47:00 

Not  good  1 

14 

24 

V.A.D 

1 5:49:00 

15:5  1:00 

.-Mt;  35,  41,  50.  60,  70,  100  1 
200.  300  at  1 rev/alt. 

14 

25 

747 

1 5:52:00 

1 5:54:00 

Azimuth  wrong,  bad  run  ! 

14 

26 

747 

15:57:00 

15:59: 00 

350  at  8 

Good  1 Best  today  | 

14 

27 

727 

1 6:02:00 

1 6:03:00 

Not  very  good  ! 

14 

28 

707 

16:05:00 

16:06:00 

340  at  8 

Looks  better 

1 ‘‘i 

29 

UC-9 

1 6:06:00 

16:07:00 

14 

30 

340  at  8 

Not  a run 

14 

31 

DC-8 

1 6:20:00 

16:23:00 

14 

32 

DC -9 

1 6:23:00 

16:24:00 
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Appendix  B 

COMPUTER  PROGRAM  DOCUMENTATION 

This  appendix  describes  the  computer  programs  developed  and  utilized 
during  the  research  program.  The  discussion  of  the  off-line  computer  soft- 
ware package  includes  flow  charts  of  the  basic  routines  and  subroutines. 


ROUTINE  ARIASK 


Finds  position  and  velo< 
of  points  in  each  normal 


Compute  nud  time  for  hi^h  velocity  points  of  scan 

■ ] 


Print  Headini; 

■ i 


r 


ROUTINE  BDOWN 


Separates  16  bit 
words  written  by 
SEL  computer 


mp 


computcii  circulation. 


^ ENTERS 
CENTER 


Are 

the  number'v 
y/  of  hijih  velocity 
/ points  leas  than  ^ 
^he  points  neressary  t< 
\ determine  a / 
vortex  / 


For  aortinjt  high  velocity 
points  as  a function  of  speed 


For  sorting  high  veloiity 
points  as  a function  of  intonsitv 


or  IS  the  number  of 
•bigh  velocity  points 
^ les  s than  S 


Set  up  Scatter  plot  with  symbols  repr 
senting  tntensilv  or  velocttv  magnitui 
depending  on  vortex  determination  ci 


r 


no  thrtsigh  high  velocity  point' 


/ this  point 

the  anale  aectorN 

to  be  excluded  / 

'y'  No 

Compute  circulation  based  on 
U8ini(  distance  from  centroid 
from  laser  through  point  and 
sight  apeed  at  point 

thi  s point 
to  line 
line  of 

__ 

Save  maximum  c irculaiion 
Call  it  (GAMA) 

Compute  circulation  baaed  on  thia  point 
uain){  diatanc  e from  centroid  to  point 
and  computed  apeed  at  point  aasummK 
circular  rotation  abmit  centroid 


Sum  circulation  | 

the  average  will  be  called  (GAMB)  j 

400  END  DO 


t ompute  a 
(GAMB 


ROUTINE  COEFF 


Computes  Fourier  coefficients 
of  VAD  sijjnal 


Fourier  coefficient  for  N harmonics 
where 

(N-1)  is  harmonic  number 
P is  number  of  points 
Y.  is  speed  at  point 
a.  is  azimuth  for  point 


t 


WINDTR 


ROUTINE  HANDLE 


Mam  logic  (or  program. 
Organises  data  to  be 
processed  by  other  nx>r 


/ 


Write  data  ironj  card  one 


L 


Return  to 
routine  ^ 
'v  inputs' 


InitMhre  record  pi'biti.in  . '‘UMer  rtnd  tin.f  « ounter 


i TAREA1>^> 


t 


•SLL  ‘ omputer  lape 


< C.,11  BDOWN  I'-P''  SEL  I 

y / i words  \ 


\ 


i Read  Card  Two  , 


( MRipole^  Hnic,  for  start  of  averaging  in  serond.<», 
lr«»m  hours,  minutes  and  seconds  I 


( oriputes  time  ol  day.  in  sec  ond.s,  from 
hours,  minutes  and  seconds 


Changes  start  and  end  frequency  for 
Spec  t rum  analvxer  from  mega  Mr  to  Mr, 


I-ooks  f -r  air.  raff  type  tn  table  to  get  wing  span 


t on»pute  s c orrelation  radius  from  wing  span  j 

(This  IS  used  in  the  vortex  detection  algorithm.)  i 


h<‘  accrplahli  putnln 


N 


I hr  nunibr  N • 

• plahU 
Ir  r iSa  n * hr 
if»'r  n^i  I -n 


I 

' r r r <r  n.r 


S*  w • ••  r pmni  ihr 

trrqurn«v  Bir  Number, 
lrtrn«iev  in  Frrquertiy  Bi 
Ran^e  »n  Focal  P<mm 
Ki  an  Angle  an 
of  Pf-int 


C»o  to  2 ^0 


Set  altitude  flag  (LDALT)  to  one  (change  in  altitude) 


Set  flag  uRROR)  to  indicate  there  ia  not  too  many  conical  acanpointi 


ave  for  this  altitude 
horizontal  wind  speed  and 
direction 

vertical  wind  speed 


U 


Go  to  350 


r 


routine  plottr 


Plots  vortex  tratks 
on  SC4020  plotter 


^ Enter  ~J 
PLOTTR 

\ / 


Y 


^Print  title,  run  number  and  time 


|~ Pact  arrays  with  positions  of  starboard., 
I port,  and  unknown  type  vortices 


For  writinK  tracks  I 
•>n  mapnetic  tape  , 


print  vortex  track  positions  \ 


\ 


Are 

■ / SC  41^20  plot^ 

/ or  printer  plotr^ 
to  be 

made  /' 

\ ,/ 

■y 

Yes 


Return  to 
HAND  I K .. 


A re 

there  SC.4020  plot). 

\ reCuested 


Are 

• hr  stairs  to  b^^i*. 
. coTnpuled  y 
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AD 

A0443I8 

a 

1 

\ 

1 

1 

1 

1 

\ 

1 

ROUTINE  PROFIL 


Project  wind  velocity  vectors  on  to  vertical 
cross  runway  and  down  runway  planes 


t 


Plots  wind  profile  as  a 
function  of  height  on  the 
SC 4020  plotter 


Return 


ROUTINE  SFLIP  Flips  sign  of  speed 

of  VAD  to  simulate 
tran  slator 


-36 


ROUTINE  SHEAR 


Computes  and  prints 
wind  shear 


■ DO  through  wind  vectors 


Compute  SHEAR  as  a function  of  altitude  in 
Wind  speed 

Wind  speed  of  vertical  cross  runway  plane 
Wind  speed  of  vertical  down  runway  plane 


These  are  computed  based  on  average  shear  between 
adjacent  pairs  of  wind  vectors.  The  shear  is  also 
calculated  from  the  derivative  of  the  power  law  curve 


-3 


ROUTINE  STATS 
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Print  the  accumulated  means 


End  DO 


Print  header  for  standard 
deviations 


DO  through  the  number  of 
minutes  averaged 


Print  standard  deviation 
for  one  minute  averaging 
periods 


End  DO 


DO  through  the  number  of  ( 37Q  ] 
minutes  averaged  V j 


Print  standard  deviations 
for  accumulated  averaging 
periods 


J (37^  End  DO 

# 

/ mean s and  \ 

1 standard  ' 

1 deviations 

\ on  unit  10  ^ 

B-4^ 


ROUTINE  WINDTR 


Computet  wind  from  VAD 
uilng  Fourier  coelficientt 


^ Enter  7 
WINDTR 


is  \ No 

Fourier  algorithm  • 


f ind  cut  off  velocity  and  count  number 
of  zero  (nonacreptable  points) 


i 


Compute  Zero  factor  for  correcting 
hnririinlrtl  velocity 


I 


Call 

COEFF 


For  computing 
Fourier  coefficients 


ents  ^ 


ompute  vc-ind  direction 


Compute  horizontal  wind  speed  with 
/•TO  correction 


Compute  raw  vertical  wind  speed 


Camiputc  vertical  wind  speed  with  a 
c or  r ec  tion  ior  the  unacc  op  table  data 


rciinp'ile  pcTientage  of  second  aiid  third 
harmonic  s and  phase  angles  I 


Print  results  of  the  Fourier  a 


Igorithm  ^ 


/ Is  si 

/wave  algorithm^ 
to  he  used 


\ 


For  computing  winds 
with  a least  squares 


ROUTINE  WPROFL 


Comp'ites  wind  profile 
as  a function  of  height 


\ Enter  / 
WPROFL 


Fit  a power  law  curve  to  horizontal  wind  velocity  at 
varying  heights  using  a least  squares  technique 


Fit  polynomial  curve  to  horizontal  wind  direction 
at  var^dng  heights  using  a least  squares  technique 


Print  coefficients  for  the  power  law  curve 
and  the  polynomial 


ROUTINE  WRT21 


Writes  vortex  tracks 
on  unit  21  for  use  by 
Data  Management  program 


WRT21 


Return 
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Appendix  B (Continued) 

In  addition  to  the  flow  charts  presented  on  the  previous  pages,  the 
following  library  routines  are  used  in  the  VAD  Vortex  Track  Program; 
FPLOT,  SOLVIT  and  LSQPF.  Rather  than  present  detailed  flow  charts  for 
these  general  library  routines,  they  are  described  below  in  terms  of  the 
purpose,  function  and  usage, 

• Subroutine  FPLOT 
Ptirpose 

FPLOT  collects  and  stores  data  as  they  become  available,  and  upon 
signal,  produces  a printer  plot  in  practically  any  orientation  and  size. 
FPLOT  should  be  regarded  as  a general  purpose  output  routine  for  dis- 
playing output  data  in  graphical  form. 

Usage 

CALL  FPLOT(Ml.  IPNT,  AR,  LR,  ISTOP,  NCU.  NCMAX.  VI,  V2) 


Description  of  Parameters 

Ml  is  the  size  of  the  main  storage  array,  and  it  should 

never  be  larger  than  800,  which  corresponds  to  400 
points  to  be  plotted. 

IPNT  is  a counter  initialized -usually  IPNT  = 0 - in  the 

calling  program.  It  is  incremented  by  2 each  time 
a new  data  point  is  entered  in  AR. 

AR  is  the  main  storage  array.  It  should  be  in  a dimension 

statement  in  the  calling  program.  For  example, 
DIMENSION  AR(800). 

LR  is  an  array  of  bytes  used  to  hold  the  curve  number. 

It  should  be  dimensioned  for  Ml/2.  The  t^ipe  declara- 
tion LOGICAL*!  LR(400)  should  be  in  the  calling  program. 

D-56 


r 


ISTOP 


NCU 

NCMAX 

VI 

V2 

IT(n 

1T(2) 

IT(3) 

IT(4) 


is  the  flag  used  to  signal  that  all  data  has  beep 
entered.  ISTOP  = 0 causes  data  to  be  stored.  If 
lSTOP  = -l  and  NCU=NCMAX  the  program  immedi- 
ately branches  to  the  plotting  section.  If  ISTOP  = l 
and  NCU  = NCMAX  a data  point  is  stored  and  then 
the  plotting  section  is  entered. 

is  the  curve  number  for  which  data  are  being  entered. 

It  must  be  a positive  integer  less  than  21. 

is  the  number  of  curves  to  appear  on  the  graph.  This 
is  the  largest  value  NCU  will  have. 

is  the  horizontal  coordinate  of  the  data  point  to  be 
plotted. 

is  the  vertical  coordinate  of  the  data  point  to  be 
plotted. 

causes  the  standard  horizontal  size  of  120  positions 
to  be  used.  A positive  integer  will  cause  the  horizontal 
size  to  be  modified.  Ma.ximum  horizontal  size  i.s  120 
positions. 

A zero  causes  the  standard  vertical  size  of  58  positions 
to  be  used.  A positive  integer  will  cause  the  vertical 
size  to  be  modified.  Maximum  vortical  size  is  238 
position  s. 

A zero  causes  the  standard  horizontal  positive  direction 
to  the  right  — to  be  used.  A 1 causes  reversal. 

A zero  causes  the  standard  vertical  positive  direction  — 
up  — to  be  used.  A 1 causes  reversal. 


LC(1  to  201  Any  symbols  placed  here  will  be  used  as  plotting 
cha  racters. 


R ema  rk  s 


As  many  as  20  different  curves  may  appear  on  a single  graph. 

A maximum  of  400  points  can  be  plotted  on  one  graph. 

The  number  of  points  in  each  curve  is  arbitrary,  but  the  total  nvmiber 
of  points  must  not  be  exceedetl. 

The  order  in  which  the  points  of  different  <urves  is  eomputed  is 
a rbit  ra  ry. 

I'ln;  graph  size  is  adjustable  up  to  120  x 2 58  positions. 

Tlie  standard  size  i s 1 20  x “58  positions  — one  page. 

The  stanfiard  liorizonta!  anfl  vertii  al  positive  directions  can  be 
independently  i hanged. 
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* Scaling  o£  data  is  automatic,  the  best  even  scale  factor  bein(j  selected 

' to  maximize  iJraph  resolution  within  the  available  space. 

i 

i A border  of  XXXXXXX  is  automatically  supplied. 

• The  background  grid  is  arranged  so  that  one  of  the  major  divisions 

1 will  pass  through  zero.  The  background  grid  is  marked  only  at  the 

j intersections  of  the  major  divisions 

i If  the  coordinate  of  zero  appears  on  a graph,  a solid  line  appears  so 

j that  this  major  division  stands  out. 

' Numerical  values  are  supplied  for  all  major  divisions. 

I The  program  plots  vertically  or  horizontally  with  equal  ease. 

When  collecting  data  for  several  curves  to  appear  on  a single  graph, 
no  plots  are  made  until  the  signal  is  given  and  the  last  data  point  for 
the  last  curve  has  been  stored. 

[ Should  one  inadvertently  try  to  store  more  data  than  storage  has  been 

[ set  aside,  data  storage  is  bypassed.  Only  the  data  stored  are  used  in 

making  the  graph. 

' A curve  with  one  point,  or  a curve  of  a constant  can  be  plotted, 

^ Stored  data  are  undisturbed  and  available  for  later  use. 

i 

f 

• Subroutine  SOLVIT 
P u rp  o s e 

Compute  a highly  accurate  solution  to  the  matrix  equation  AX  = B 
where  A is  an  N x N matri.x,  and  X and  B are  vectors. 

Usage 

CALL  SOLVIT  (A,  N,  M,  B,  X,  ACC.  MAXIT,  IT,  IN,  W) 


Description  of  Parameters 


A 

N -X  N matrix  of 

N 

numbe  r 

of  rows 

M 

numbe  r 

of  rows 

B 

right -ha 

nd  side 

X 

sol  ution 

vecto  r 

coefficients 
of  A 

of  array  in  which  A is  stored 
of  equation  AX  = B 
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ACC 

accuracy  desired 

MAXIT 

maximum  number  of  iterations  to  be  tried, 
be  sufficient.)  If  MAXIT  = 0,  no  iterations 
performed. 

(10  should 
will  be 

IT 

internal  indicator  set  by  program  to; 
= - 1 if  matrix  is  singular 

= 0 if  convergence  failed  within  MAXIT  iterations 
= k the  number  of  iterations  it  took  to  converge 

IN  to  be  set  = 1 for  first  entry  to  the  subroutine 

= 2 for  subsequent  entries  if  A has  not  been 
changed. 

W Nx(N+5)arrayusedas  working  storage. 

Method 

The  Gaussian  elimination  method  with  partial  pivoting  is  used  to  de- 
compose the  matrix  to  upper  and  lower  triangular  matrices.  (These  tri- 
angular matrices  are  saved.)  From  the  residual  vector  B - AX.  a correction 
vector  is  computed  and  added  to  X.  The  process  is  repeated  until  the  re- 
sidual vector  is  within  the  tolerance. 


The  initial  scaling  of  each  row,  to  bring  the  maximum  element  between 
1/2  and  1.  is  done  using  factors  equal  to  integral  powers  of  2 to  eliminate 
round  off  errors  which  otherwise  might  have  been  introduced. 

The  method  is  due  to  Clcve  Moler.  ,TPL,  presented  at  SHARE  XXV. 
Chicago. 

Storage  Requirements 
1 1 6 location s. 


R est  rictions 


In  most  cases,  eight  digits  accuracy  can  be  attained. 
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W must  not  be  used  between  calls  to  SOLVIT  when  solvin'^  for  multiple 
right-hand  sides.  Matrix  A is  destroyed. 

Note 

SOLVIT  is  only  a driver  for  the  deck  SIME  (entry  points  SIMEQA  and 
SIMEQB);  thus  two  decks  comprise  this  subroutine.  In  addition,  this  sub- 
routine uses  the  library  subroutine  LOG2. 

• Subroutine  LSQPF,  LSQPFl 
Purpose 

This  program  will  fit  polynomials  of  order  one  through  seven  to  n 

given  points  (Xj  , Y),  (X2,  Y^), ^^n’  ^n'  method  of  least  squares. 

The  order  and  spacing  of  the  points  is  immaterial  and  the  points  need  not  all 
be  distinct. 


U sage 

a.  CALL  LSQPF  (X,  Y,  K,  N,  M,  C.  lERR)  or 

b.  CALL  LSQPFl  (X.  Y,  K,  N,  M,  C,  lERR) 


Description  of  Parameters 

X for  both  entries  X is  the  first  location  of  the  block 

of  X^  values  (real  array! 

Y for  both  entries  Y is  the  first  location  of  the  bb  ck 

of  Y.  values  (real  array) 

K for  both  entries  K is  zero  if  residuals  are  not  needed; 

is  nonzero  if  residuals  are  wanted  (integer) 

N for  both  entric^s  N is  the  number  of  points  to  be  fitted 

(ini  ege  r) 

M (or  entry  (a),  M is  the  highest  order  poUnoniiai  desired. 

The  routine  will  compute  polvnomial  fils  of  degrees  one 
through  M and  return  all  coefficients.  For  entry  (!>), 
the  routine  computes  a pol^-nomial  fit  of  degree  M only. 

2 < M ■ 7 (integer) 
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f 


c 


lERR 


for  both  entries,  C is  the  first  location  of  the  block 
where  the  coefficients,  standard  deviation,  and  residuals 
will  be  stored. 


for  entry  (a),  the  dimension  of  C must  be  M(M+5)/2  if 
K = 0,  and  M(M+5)/2  + N * M if  K / 0.  The  coeffi- 
cients, a^,  the  standard  deviation,  a,  and  the  residuals, 

R are  stored  in  the  following  order  beeinning  at  loca- 
J' 

tion  C : 


A<*>, 

o 

A</>, 

a<2), 

O 

A«>, 

Af. 

— A<‘^>, 
m 

Rj'  \ 

r<^»,  -- 

n 

The  location  of  of  polynomial  degree  P(1  ^ P < M) 

is  C(I)  where  I = P(P+5)/2  - fP+1).  The  location  of 
Rj  of  polynomial  degree  P(1  P ^ M)  is  C(I)  where 

1 = M(iM+5)/2  + N(P-1)+1. 

For  entry  (b),  the  dimension  of  C must  be  (M+2)  if 
K = 0.  and  (M+2+N)  if  K ^ 0.  The  coefficients  are 
stored  in  C(l)  through  C(M+1),  the  standard  deviation 
in  C(M+2),  and  the  residuals  in  C(M+3)  through 
C(M+2+N).  (real  arrayl 


for  both  entries,  lERR  is  the  error  flag;  it  will  contain 
0 if  no  error  occurred,  and  some  number  if  an  error 
ha.s  occurred.  There  are  two  reasons  for  an  error; 
one  of  the  pivotal  elements  is  zero,  or  an  overflow,  has 
occurred.  (integer) 


Methorl 

The  method  of  solution  is  as  follows; 

Given  a set  of  n points.  (X..  Y.).  we  require  sets  of  coefficients,  a., 
such  that 


I 


r.  -h  1 


r 


is  the  best  least  squares  fit  to  over  all  i.  This  is  done  by  first  normal- 
izing the  X.  to  the  interval  (-1,+1)  by  use  of  the  formula 


max  1 

' |X  1,  . . 

..IX  1 \ 

1 

J o 1 

1 n 

We  can  then  work  in  terms  of  Chebychef  polynomials,  T,  , where 


T =1 

0 


Ti  = X- 

Tk  = 2 X. 


for  k = 2(1  )p.  We  now  seek  h such  that 


E ‘j  b (») 

j=0 

gives  the  best  least  squares  fit.  The  residuals,  r!',  are 

k 

= 23  [®j  <^i>'’]  ■ '^i' 

j=0 

The  unbiased  estimate  of  the  standard  deviation  of  the  k^''  order  polynomial 
fit  is 

a(k'  = 


R t'st  rictions 

The  valut^of  M must  be  in  the  range  2 through  7.  Care  should  be 
exercised  in  using  this  routine  since  considerable  loss  in  accuracy  of  the 
l oefficients  is  sometimes  involved. 
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Appendix  D 

REPORT  OF  ENT  IONS 

Under  this  contract  the  Lockheed  LDV  system  was  modified,  calibrated 
and  tested  in  an  airport  environment.  A base  of  wind  and  wake  vortex  meas- 
urements were  obtained.  The  contract  objectives  were  met  and  no  inno- 
vation, discovery  nor  invention  was  made. 
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